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Abstract 

Cerenkov images of air showers can also be classified using multifractal and wavelet 
parameters, as compared to the conventional Hillas image parameters. This new 
technique was applied to the images recorded by the cameras of the stereoscopic 
imaging air Cerenkov-telescopes operated by the HEGRA collaboration. With re- 
spect to the identification of particles, the performance of multifractal and wavelet 
parameters was examined using a data sample from the observation of the active 
galaxy Mkn 501 that showed a high 7-ray flux. The multifractal parameters were 
also combined with the Hillas parameters using a neural network approach in order 
to further improve the 7/hadron-separation. 

Key words: imaging air Cerenkov-technique, 7/hadron-separation, multifractal 
parameters, wavelets, neural networks 
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1 Introduction 

Over the last decade, imaging air Cerenkov-telescopes have emerged as the 
most powerful tools for 7-ray astronomy in the TeV energy range. The imaging 
of air showers allows a separation of 7-ray induced air showers and hadron- 
induced air showers, and therefore a strong reduction of the cosmic-ray back- 
ground ([1,2]). The performance of such instruments can be boosted further 
by combining several Cerenkov-telescopes to form a stereoscopic system, pro- 
viding several simultaneous views of an air shower. For example, the HEGRA 
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system of Cerenkov-telescopes ([3-5]) achieves an angular resolution of 0.1° 
for individual photons, and an energy resolution of better than 20%. 

In the traditional approach, the distribution of the Cerenkov light in the cam- 
era is characterized by an ellipse, whose semi-axes are derived by the calcu- 
lation of the second moments of the amplitude distribution after a suitable 
procedure to clean the image of singular high amplitudes at a large distance 
from the actual projection of the shower ([6,7]). The width of this ellipse very 
efficiently distinguishes between electromagnetic and hadronic showers. To re- 
move the dependence of the image width on zenith angle, total amplitude and 
distance of the shower core to the telescope, one can use the scaled width ([8]) 
where the width parameter is divided by the Monte-Carlo expectation value for 
7-rays. For a system of Cerenkov-telescopes, the scaled width is averaged over 
all triggered telescopes. The resulting mean scaled width of the image ellipses 
peaks at a value of 1.0 for 7-originated showers and shows excellent perfor- 
mance with respect to separating 7-rays from the hadronic background ([8]). 
Thus, cuts on mean scaled width enhance the significance of a 7-signal strongly. 



In an alternative approach, one may visualize the cascading process in which 
a shower develops as a series of self-similar processes and hence one expects 
that the distribution of Cerenkov light exhibits self-similar structures, as was 
pointed out in [9] . The possibility exists that in a shower of reasonable size the 
superposition of the Cerenkov light from all sub-showers results in a smooth 
image in which the these structure has been averaged out, especially when 
taking account of the fact that the Cerenkov light is subjected to scattering 
processes. Therefore, the Cerenkov light distribution recorded by the camera 
does not necessarily have to have fractal properties. 

Mult if r act al analysis aims at describing this self-similarity quantitatively. The 
scope of the article is to give a brief introduction to multifractal and wavelet 
analysis (section 2), to describe how these methods were applied to the images 
recorded by the HEGRA-cameras, and to investigate their properties in deter- 
mining the nature of the primary particle by comparing their performance on a 
Mkn 501 data sample with the approach by Hillas (section 3). An attempt was 
made to combine both Hillas parameters and multifractal parameters using 
neural networks (section 4). 



2 Multifractal Analysis 
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2. 1 Multifractals 



The quantities derived in multifractal analysis describe the behaviour of pat- 
terns with respect to scale transformations. We assume that the pattern - 
i.e., the Cerenkov image - is described as a density A(xl) defined on a two- 
dimensional grid, corresponding to the amplitudes measured in the photomul- 
tiplier pixels of a camera. In a first step, the distribution is normalized such 
that Yj{xi} A(xt) = 1. Starting point of the analysis is the calculation of a 
multifractal moment Z q (£) of order q, defined on a length scale £ ([10]) 

z q (£) = '£p(x,e) q with p(x,£)= (!) 

{x} g box of size i 

centered at x 

Here, the image is divided up into a set of boxes of characteristic scale l\ the 
non-overlapping boxes have identical shapes and cover the image completely. 
The amplitudes A(xf) within a box of size £ at a lattice site x are summed 
up, yielding p(x,£). These so-called box-amplitudes are exponentiated with a 
number q, that assumes integer values, and successively the summation over 
a set boxes is formed. 

With varying lengthscale £ one expects for a genuine fractal distribution a 
behaviour like 

Z q {£)~£ T ^\ (2) 

from which the first multifractal parameter, the mass exponent r(q), may 
be derived as the slope of the line through the points (Z q (£)\£) in a double- 
logarithmic plot. The order q of the moment defines how structures that are 
less bright and prominent contribute to the value of Z q (£). The singularity 
exponent a(q), defined as 

describes how r(q) changes with changing order q. 
2.2 Wavelets 

Wavelet analysis is closely related to multifractal analysis: it describes the 
scaling property of differences in amplitude in adjacent boxes rather than the 
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amplitudes themselves, i.e. the derivative of the amplitude distribution on dif- 
ferent baselines. Alternatively one may think of wavelet analysis as a Fourier 
expansion with a discrete set of basis functions that are generated by scaling 
from a so-called mother wavelet ([11]). 

In analogy with multifractal analysis ([12]), one defines a wavelet moment 
whose behaviour under scale transformation is determined by an exponent (3. 
For a one-dimensional density, the wavelet moment is given by 

W q (£) = £ \ P (x,£) -p(x + £- e x ,£)\ q W q (£) ~ £^\ (4) 

{x} 

The relevant quantity for the derivation of a wavelet moment is the differene 
in integrated amplitudes between one box of linear extend £ and the neigh- 
boring box in the direction indicated by the unit vector e x . This corresponds 
to the convolution of the amplitude distribution with a wavelet (in this case 
the simplest so-called Haar wavelet), anlyzing the distribution with at spatial 
resolution I. The sign of the wavelet function is indicated pictorially by (H — ). 

In order to characterize a two-dimensional distribution on a scale £ by its 
wavelet moments, a set of three base wavelets , 1) , +) is required 
([12]), each yielding a wavelet moment and a corresponding exponent. 

2. 3 Application to the Shower Images 

The HEGRA-Cerenkov-telescopes ([13]) are equipped with cameras ([14]) con- 
sisting of 271 pixels of 0.25° diameter arranged on a hexagonal lattice. The 
telescopes are triggered by a coincidence of two pixels with at least 6 to 8 
photoelectrons each. Images contain about 100 photoelectrons per TeV of 7- 
ray energy. Therefore, images near the trigger threshold (0.5 TeV) contain 
about 50 photoelectrons, typical images about 100 photoelectrons, and the 
largest gamma-ray images about 2000-3000 photoelectrons. Night-sky back- 
ground light as well as electronics noise and ADC quantization cause a typical 
noise of 1 photoelectron rms in each pixel. For the normal Hillas-type analysis, 
only pixels with signals of at least than 6 photoelectrons are used, or pixels 
with at least 3 photoelectrons provided that they are adjacent to a pixel with 
at least 6 photoelectrons. 

The hexagonal lattice of the cameras complicates the fractal analysis slightly. 
Unlike in the case of a square lattice ([9]), it is not advisable to form "boxes" 
of 1,4 = 2x2, 9 = 3x3,... pixels, since the resulting "boxes" have a preferred 
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axis. Instead, hexagonal cells consisting of 1, 7, 19 and 37 pixels were defined; 
the size of one cell is described by the extent of the hexagon along one side, 
thus £ assumes values of £ = 1, 2, 3 and 4, as it is shown in Figure 1. 




Fig. 1. Definition of the cells of pixels for £ = 2 (top row), £ = 3 (middle row) 
and £ = 4 (bottom row), superimposed to a shower image containing about 2000 
photoelectrons (left column) and about 100 photoelectrons (right column). The cells 
corresponding to £ = 1 are the individual pixels. The ampitude in a given pixel is 
indicated by the block area. The thin line indicates the border of the camera. 

Another complication arises from the noise contributions. There is a finite 
probability to have a very small signal in a pixel, which would make moments 
with negative index q diverge. Since the amplitude distribution is more or 
less continuous - the peaks corresponding to 0, 1, 2, ... photoelectrons are not 



5 



resolved individually - such pixels could only be removed by a more or less 
arbitrary cut. Instead, it was decided to limit the analysis to positive values 
of q, in the range from 2 to 8. The values q = and q — 1 yield trivial results; 
q = simply counts the numbers of cells and the moment for q — 1 simply 
reflects the normalization of the amplitudes. 



3 Results 

The analysis methods described above were applied to a subset of the data 
resulting from the observation of the active galaxy Mkn 501 in a high state 
of activity in the course of the year 1997 ([4], [5]). The data were taken at 
zenith angles between 10° and 26°. During these observations, the source was 
displaced by 0.5° in declination with respect to the optical axis of the tele- 
scopes; a region displaced by the same amount in the opposite region served 
to evaluate backgrounds under the signal. The sign of the displacement varied 
between 20-minute runs, in order minimize systematic effects. The distance of 
1° between the signal region and the background region is large enough that 
events reconstructed stereoscopically can be assigned unambiguosly. 

After a preselection keeping only events with shower axes within 0.5° from 
the source or from the center of the background region, the data set contained 
approximately 29.000 events originating from the region centered on the source 
position opposed to 25.000 hadronic background events. Images showing a 
total amplitude of more than 40 photoelectrons were included in the analysis. 

3. 1 Distributions 

Figure 2 shows the distributions of the multifractal parameters r 2 and a 2 for 
7-rays as well as for cosmic rays. The distribution is very well described by a 
Gaussian function. The multifractal parameter r indicates how fast the num- 
ber of cells of pixels that carry a significant amount of the total amplitude 
increases, if the size of the cell is decreased. Images originating from a 7-ray 
yield smaller values for r, because these images are well contained in a narrow 
region of the camera in contrast to the fuzzy images produced by hadronic 
showers. For these hadronic images, one observes a fast increase in the num- 
ber of cells with decreasing cell size, and therefore the r-distribution peaks at 
larger values. 

The parameter a indicates how fast r changes with increasing order q. Consid- 
ering the images of a hadronic shower, they are spread out over a large region 
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Fig. 2. The distributions of the mass exponent r 2 (left) and the singularity exponent 
a 2 (right) for 7-rays (open circles) und hadrons (filled circles). 

of the camera and thus the different cells of pixels carry only small amplitudes, 
i.e. numbers close to zero in contrast to a 7-induced shower, which generates 
amplitudes that are relatively large in the various cells of pixels. With increas- 
ing q the moments Z q (£) decrease rapidly in the case of images produced by 
hadrons, resulting in a fast increase in the slope r(q) with growing order q. 

The distribution of the wavelet parameter /3 2 is shown in Figure 3. Images of 
hadronic origin show larger fluctuations on the short length scales compared 
to an image produced by 7-rays, that may in fact be described by a smooth 
gaussian distribution. This results in an increase of the wavelet moments W q (£) 
with decreasing length scale. Therefore, the parameter (5 assumes larger values 
for hadronic events compared to 7-rays. 



3.2 Separation Properties Compared with the Hillas-P ammeters 

In order to describe the effectiveness of the different parameters with respect 
to distinguishing between 7-rays and the hadronic background, cuts on these 
parameters were applied. From that, efficiencies k and Q- values were derived: 

K 7 = Nf^J/Ny, KHadron = ^Hadron/^Hadron and Q = K y / ^K^dron- The num- 
ber iV 7 was determined by statistically subtracting the numbers of the events 
in the source region and in the background region. The Q-value describes by 
how much the significance of a weak signal, calculated as in [15], is enhanced 
by the cuts. The analysis was first applied to all telescopes individually. 
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wavelet exponent /? 2 

Fig. 3. The distribution of the wavelet parameter 02 for 7-rays (open circles) and 
hadrons (filled circles). 

Table 1 lists mean values and standard deviations of the distribution for con- 
ventional and multifractal image parameters for the 7-rays as well as for the 
hadron-distributions, the 7-efnciencies re 7 , hadron-emciencies /tHadron and Q- 
values for a cut that yields the largest Q- value in descending order with respect 
to the Q- value they archieve. As already evident from Figure 2, multifractal 
parameters can be employed to enhance the significance of 7-ray signals, but 
seem not to be competitive with the conventional parameter scaled width. 





scaled 
width 


T2 


012 


scaled 
length 


ih 


(XT) 


1.04 


1.18 


1.12 


0.99 


1.45 


°l 


0.18 


0.24 


0.24 


1.26 


0.50 


^Hadron^ 


1.58 


1.47 


1.40 


0.19 


1.81 


Hadron 
X 


0.38 


0.24 


0.26 


0.29 


0.56 


cut(°^ 


sw < 1.1 


r 2 < 1.25 


a 2 < 1.25 


si < 1.2 


f3 2 < 1.9 




0.663 ± 0.044 


0.722 ± 0.050 


0.797 ±0.056 


0.710 ± 0.006 


0.879 ± 0.067 


K Hadron 


0.113 ±0.004 


0.241 ± 0.005 


0.338 ± 0.006 


0.404 ± 0.054 


0.635 ± 0.008 


Q-value 


1.82 ±0.12 


1.47 ±0.19 


1.37 ±0.09 


1.12 ±0.08 


1.10 ±0.08 



Table 1 



Description of the distributions of the various parameters, including efficiencies and 
the Q-value for an optimal cut, i.e. for maximum enhancement of significance. The 
values quoted refer to events detected with the telescope CT3; values resulting from 
the analysis of images in the other telescopes are consistent within statistics. 
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3.3 Dependence on q 

An interesting point is the dependence of the parameters on the index q, es- 
pecially how the separation performance of the parameters is influenced by 
different weighting exponents. Figure 4 illustrates how the parameters r and 
a change with q. Depicted are the mean values for the 7-ray and hadron dis- 
tributions a well as their widths resulting from a Gaussian fit. 




12345 6 789 1234 5 6789 

multifractal irdex q multifracta irdex q 



Fig. 4. Mean values and standard deviations of the multifractal parameters r (left) 
and a (right) as a function of the index q for 7-rays (open circles) and hadrons 
(filled circles). 

The parameters show the expected dependence on q. However, their separa- 
tion performance is virtually unaffected by a particular choice of q - in all 
cases, the separation between 7-rays and hadrons corresponds roughly to the 
width of the distributions. The reason is that parameters of different order q 
turn out to be strongly correlated. 

Furthermore, the values of the multifractal parameters are affected by the 
position of the shower image relative to the pixel cells, clearly an effect of the 
relatively large cells. The value of the moment Z q (£) depends on the position 
of the image relative cells of pixels especially for the largest cell. Also the 
points in the double-logarithmic plot of Z q {£) against £ do not lie on a line 
as expected from equation 2, but rather on a concave curve. For that reason 
the slope r of a line through these points gets steeper if one includes only the 
moments Z q (£) on the short length scales. 
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3.4 Correlations 



An interesting question is to which extent the multifractal parameters and 
the Hillas parameters are correlated. If multifractal parameters are genuinely 
sensitive to fluctations within the image, a weak correlation would be expected. 
If, on the other hand, the multifractal parameters merely provide a different 
parametrization of the global shape of the image, one would expect strong 
correlations. 

In Figure 5 the correlation between the multifractal parameter r and the 
Hillas-parameters scaled width and concentration is depicted for the hadronic 
background. The mean value and the standard deviation of the Hillas param- 
eter is plotted as a function of the value of a multifractal parameter. 




mass exporert t 3 mass exporert r 



Fig. 5. Correlation of the multifractal parameter r with the Hillas-parameters scaled 
width (left) and concentration (right). 

Clearly, the parameters are not independent of each other. The correlation 
of the multifractal parameter r with scaled width directly follows from the 
arguments outlined in section 3.1. Furthermore, the correlation of r with con- 
centration shows an accumulation of 7-candidates at small values for r and 
large values for concentration, whereas hadron-candidates amass at large val- 
ues for r and small values for concentration. 

3.5 Averaging over Multiple Telescopes 

In the conventional analysis using the (scaled) Hillas parameters, separation 
properties are enhanced by averaging over the images taken from the same 
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shower by the different telescopes. The same is true for the multifractal pa- 
rameters. Averaging values of the multifractal parameters r or a over all 
telescopes that have been triggered by a given event results in a significant 
improvement of its separation properties, as can be seen in Table 2. Still, 
however, are the multifractal parameters inferior in their performance. 





cut 




^Hadron 


Q-value 


T2 


r 2 < 1.17 
fa < 1-60 


0.555 ± 0.035 
0.650 ± 0.046 


0.076 ± 0.003 
0.318 ±0.005 


2.02 ±0.13 
1 15 ±0.08 


mean scaled width 
mean scaled length 


msw < 1.10 
msl< 1.10 


0.718 ± 0.043 
0.724 ± 0.039 


0.074 ± 0.003 
0.326 ± 0.005 


2.64 ±0.16 
1.27 ±0.09 



Table 2 



Efficiencies and Q-values for an optimal cut for the averaged parameters. 



4 Neural Networks 

In order to investigate how the 7/hadron-separation may be improved in spite 
of the close correlation of the parameters, a neural network was employed. A 
detailed introduction to the functionality and theory of neural nets may be 
found in [16]. The analysis was applied to data from individual telescopes. 
Neural networks of the type 'feed-forward perceptron' having 

(1) the Hillas-parameters width, length, concentration and log(sue); 

(2) multifractal and wavelet parameter T2, T4, t$, fa, Pa, fis and log (sue); 

(3) all of the above parameters 

as input values were implemented and their classification performance was 
examined. The adjustment of the net's weights and thresholds was done using 
a data sample consisting of 1200 events with a reconstructed shower direction 
within a circle of 0.2° around the source-position as 7-rays and 1200 events of 
background data as hadrons. The events were requested to have an image-size 
of at least 40 photoelectrons. The network was adjusted in such a way that 
it produced an output value of z = for the 7-candidates and z — 1 for the 
hadron-candidates, respectively. 

The neural nets were implemented using the Je£/Vei-programme ([17]), that 
provides a number of interesting features for optimizing the adjusting process 
such as adding momentum or introducing artificial noise, in order to prevent 
the system from being trapped in a local minimum. The topology of the net- 
work was chosen in such a way that the networks consisted of three layers in 
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all three cases, but the number of nodes within a layer was altered until the 
net yielded the best result on the training data. 

Figure 6 illustrates the distribution of the output z of the neural network, 
clearly showing an excess of events with low values for z corresponding to 
7-rays. 




a> 0.2 0.4 0.6 0.8 

neural network output z 



Fig. 6. The distribution of the output z of the neural net using Hillas- as well as 
multifractal parameters; filled circles correspond to the hadronic background, open 
circles to 7-rays. Both distributions have been normalized to unit area. 





cut 




^Hadron 


Q-value 


1. 


z < 0.38 


0.70 ± 0.05 


0.15 ±0.004 


1.79 ±0.12 


2. 


z < 0.34 


0.54 ±0.04 


0.12 ±0.003 


1.59 ±0.11 


3. 


z < 0.26 


0.54 ±0.05 


0.07 ±0.003 


2.01 ±0.14 



Table 3 

Efficiencies and Q-values for optimal cuts on the neural net output z, comparing 
the various approaches. 

The network based on Hillas-parameters was able to reach the performance 
of the parameter scaled width, that attains a Q-value of 1.82 ± 0.12. The net- 
work using all available parameters slightly surpasses the performance of the 
network that relied on Hillas-parameters. The network based on multifractal 
and wavelet parameters falls behind, as can be seen from Table 3. 
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5 Concluding Remarks 

Image analysis techniques based on multifractals and wavelets were applied to 
images recorded by the HEGRA-Cerenkov-telescopes. A 7/hadron-separation 
by cuts on these parameters is feasible, yielding a performance similar to the 
Hillas-parameter concentration, but falling behind that of the parameter scaled 
width. Combining Hillas parameters and multifractal parameters using a neu- 
ral network, a slight improvement in performance could be achieved compared 
to the Hillas parameters alone. The gain is, however, limited by the strong cor- 
relation between the Hillas parameters and the multifractal parameters. 

A general - and not too surprising - conclusion is certainly that the pixel 
size of the HEGRA cameras is much too coarse to effectively characterize the 
fractal nature of showers images; rather than exploring the fine structure of 
the image, the multifractal parameters determine its overall shape, explaining 
the good correlation with the Hillas parameter concentration. 

Our overall assessment concerning the use of multifractal parameters to iden- 
tifiy the primary particle of an air shower is less optimistic than the view 
given in ([9]), in which the analysis method is applied to Monte-Carlo data 
from a simulation of the TACTIC-experiment. While the pixel size of HEGRA 
(0.25°) and TACTIC (0.31°) are similar, some basic differences between the 
two analysis should be emphasized: 

• The square-grid pixel arrangement of TACTIC simplifies the subdivision 
into pixel cells. 

• The analysis of [9] is applied to images containing at least 1800 photoelec- 
trons, whereas the HEGRA analysis was applied to images with as little 
as 40 photoelectrons. For such small numbers of photoelectrons, fluctua- 
tions in the image are dominated by photoelectron statistics, rather than 
the shower substructures. For gamma-hadron separation, algorithms which 
work only for very large images are of limited interest. The situation differs 
if one applies the technique to separate primary nuclei, which is the main 
emphasis of [9] . 

• The analysis of [9] also included multifractal moments for negative expo- 
nents q. This was possible since in the simulation Poisson-distributed back- 
ground photoelectrons were added and there is a clear distinction between 
and 1 photoelectron in a pixel. In HEGRA, with a Gaussian noise of about 
1 photoelectron rms, an ad-hoc cut would have to be used to exclude low 
amplitudes; to avoid this, only positive values of q were allowed. 

In order to really judge if one can detected the fractal structure of air show- 
ers in their Cerenkov images, one should probably employ a simulation with 
very fine pixels - such that the first set of length scales is well contained 
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within the image - and a photoelectron yield such that Poisson fluctuations 
are small compared to the shower fluctuations. The first condition implies a 
pixel size around 0.002° to 0.01°, the second an image size of 10000 to 100000 
photoelectrons. 
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